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Topological spin transport by Brownian diffusion of domain walls 
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We propose thermally-populated domain walls (DWs) in an easy-plane ferromagnetic insulator as 
robust spin carriers between two metals. The chirality of a DW, which serves as a topological charge, 
couples to the metal spin accumulation via spin-transfer torque and results in the chirality-dependent 
thermal nucleation rates of DWs at the interface. After overpopulated DWs of a particular (net) 
chirality diffuse and leave the ferromagnet at the other interface, they reemit the spin current by spin 
pumping. The conservation of the topological charge supports an algebraic decay of spin transport 
as the length of the ferromagnet increases; this is analogous to the decaying behavior of superfiuid 
spin transport but contrasts with the exponential decay of magnon spin transport. We envision that 
similar spin transport with algebraic decay may be implemented in materials with exotic spin phases 
by exploiting topological characteristics and the associated conserved quantities of their excitations. 

PACS numbers: 75.76.+j, 75.78.-n, 66.30.Lw, 75.10.Hk 


Introduction. —Spintronics, or spin-transport electron¬ 
ics, exploits spin degrees of freedom in condensed mat¬ 
ter systems to improve information processing technology 
that is traditionally based on electric charge [1]. Con¬ 
ducting materials have been used to transport spin by 
polarizing itinerant electrons, which is associated with 
undesired energy dissipation due to the electronic con¬ 
tinuum. Magnetic insulators, which are immune to Joule 
heating, provide alternative platforms to seek an efficient 
spin transport channel. Superfluid spin transport Pi¬ 
ll] has been proposed for long-ranged spin transmission 
in magnetic insulators with easy-plane anisotropy. The 
spin superfluidity, however, can be destroyed by U(l)- 
symmetry-breaking anisotropy within the easy plane. 

Topological solitons in magnetic materials are non¬ 
linear excitations that are protected by their nontriv¬ 
ial topology [6|. A domain wall (DW) in an easy-axis 
magnet is a prototypical topological soliton, which can 
store and deliver information as demonstrated in the 
racetrack memory [7|. DWs can be driven by various 
means, e.g., an external magnetic field [8], an electric 
current (in conducting systems) [9|, or heat flux [TOHlB] . 
At a hnite temperature, DWs with damped dynamics un¬ 
dergo Brownian motion due to a random force dictated 
by the fluctuation-dissipation theorem HMD; under a 
temperature gradient, Brownian motion leads to a diffu¬ 
sive transport (thermophoresis) of DWs [22] . 

In this Rapid Communication, we show that 
superfluid-like spin transport can be achieved by utilizing 
thermally-populated DWs in an easy-plane ferromagnetic 
insulator with an additional easy-axis anisotropy within 
the easy plane. Long thin ferromagnetic strips, for ex¬ 
ample, are naturally endowed with such anisotropies due 
to magnetostatic interactions See Fig. for illus¬ 

tration. A DW is characterized by its chirality q = ±1, 
associated with the sense of circulation of the magne¬ 
tization within the easy plane [6]. The chirality of a 
DW is protected in the XY ferromagnet by topology. 
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FIG. 1. (color online). An easy-a;y-plane ferromagnetic in¬ 
sulator with an additional easy-axis anisotropy in the x di¬ 
rection is sandwiched between two metals. The spin-transfer 
torque caused by the out-of-equilibrium spin accumulation 
fj, in the positive 2 direction prefers injection of DWs with 
the clockwise-rotating magnetization. The annihilation of 
these DWs generates the spin current into the right metal 
via spin pumping. In the diffusive limit of DW motion, the 
spin current decays algebraically as the ferromagnet’s length 
increases. 


and we thus refer to it as the topological “charge.” Sup¬ 
pose the ferromagnet is driven out of equilibrium by the 
spin accumulation in the positive 2 direction in the left 
metal. The induced spin-transfer torque nucleates DWs 
with the clockwise-rotating magnetization. When these 
DWs leave the ferromagnet toward the right metal, the 
magnetization at the interface rotates counterclockwise, 
which, in turn, generates the spin current into the metal 
via spin pumping. In the diffusive regime of DW mo¬ 
tion, the spin current transported by DWs decays alge¬ 
braically as in superfluid spin transport owing to the 
conservation of the topological charge. This topological 
spin transport can be inferred by measuring the drag co¬ 
efficient in a magnetoelectric circuit that was proposed 
in Ref. |S] for detecting superfluid spin transport. 

Main results .—The model system consists of a quasi 
one-dimensional easy-a:y-plane ferromagnetic insulator 
with an additional easy-x-axis anisotropy attached on 
both sides by nonmagnetic metals. In equilibrium, the 
anisotropy lays the local spin density s = sn in the xy 
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FIG. 2. (color online), (a), (b) The DWs with the topological 
charge <7=1. (c), (d) The walls with q = —1. 



plane, which allows us to parametrize its direction as 
n = (cos (j), sin (j),0). A DW is a topologically stable equi¬ 
librium texture that interpolates the two uniform ground 
states, (/) = 0 or TT. Its associated winding is characterized 
by the topological charge: 


q = 


1 

TT 


/ 


dxd^cj), 


( 1 ) 


FIG. 3. (color online). Schematic for the conservation of 
the total topological charge, (a) A pair of DWs with oppo¬ 
site charges, so that the direction of the magnetization does 
not wind the circle as shown in the right. The magnetiza¬ 
tion textnre is, therefore, topologically trivial and can be cre¬ 
ated or destroyed spontaneously, (b), (c) A pair of DWs with 
the same charge. The direction of the magnetization winds 
around the circle once, which makes the textured configura¬ 
tion topologically stable from thermal annihilation. The total 
topological charge, i.e., the net winding number, is conserved 
during interactions between DWs. 


where the integral is over the DW along the longitudi¬ 
nal X axis of the ferromagnet. Figure illustrates four 
possible DW types. 

A finite temperature causes spontaneous nucleation 
and annihilation of DWs. In the bulk, DWs are created 
and destroyed always in pairs with opposite charges as 
shown in Fig. [^a) [2l]. The topological charge density, 
p = p'^ — p~ is, thus, preserved in the bulk [Fig. [^b) 
and (c)], where p^ are the linear densities of DWs with 
q = ±1, respectively. A topological charge can be 
injected or ejected through the boundaries of the fer¬ 
romagnet. In equilibrium, the DW density is charge- 
independent; p^ ^ po (X exp(—A q/T), where Eq is the 
DW energy. 

A DW should generally behave as a particle immersed 
in a viscous medium due to its coupling to, e.g., lattice vi¬ 
brations [18] or other microscopic degrees of freedom. As 
such, it must exhibit Brownian motion at a finite temper¬ 
ature due to random forces, whose existence is dictated 
by the fluctuation-dissipation theorem m- For a con¬ 
glomerate of DWs that diffuse by Brownian motion, the 
dynamics of the topological charge density is described 
by the Fokker-Planck equation [22] : 

dtp -b 0/ = 0, I = -Dd^p , (2) 

in the absence of an external force, where I is the topo¬ 
logical charge current. In equilibrium, the density and 
the current of the topological charge are zero; p = 0 = I 
according to the reflection symmetry in the xz plane and 
the time reversal symmetry. 

The topological charge density can be injected by per¬ 
turbing the ferromagnet by the nonequilibrium z axis 
spin accumulation in the left metal, /x = /iz, which 
we assume positive, p > 0, for concreteness. The 
spin-transfer torque caused by the spin accumulation is 
T = {g'^+gLnx){fj,xn)/4Tr, where gj^ = gi-\-ig[ is the ef¬ 
fective complex spin mixing conductance associated with 


the ferromagnet/metal-f interface [IS]. The torque does 
work on the ferromagnet favoring the nucleation of DWs 
with the positive charge: W® = qgLpS/A, where q is 
the charge of the wall and S is the cross-sectional area 
of the ferromagnet. The resultant nucleation rate of the 
topological charge is F lSW/T to linear order in the bias, 
where F^ is the equilibrium-nucleation rate of DWs at 
the left interface and 5W = FF+ — W~ = gLpS/2 is the 
difference between the two works. 

The injected topological charges diffuse by Brownian 
motion and can leave the ferromagnet through the right 
boundary. The conservation of the topological charge 
leads to the steady-state current (as derived below): 

I = - -^ ^3^ 

Rl + Rr + Rb 


where 


Rl = 


2T 

f0’ 


Rr = 


2T 


Rb 


2TL 


( 4 ) 


and L is the length of the ferromagnet. We may in¬ 
terpret the topological charge current I as the applied 
“voltage” gLP (with units of J/m^) divided by the total 
“resistance” Rl + Rr + Rb (with units of J • s/m^) of 
the series circuit, which is made of the interface resis¬ 
tances, Rl and Rr, and the bulk resistance Rb- Note 
that the bulk resistance Rb is proportional to the ratio 
of the length to the cross-sectional area, L/S, which is 
analogous to the electrical resistance. 

The dynamics of the local spin density at the bound¬ 
aries injects spin current into the metals via spin pump¬ 
ing, which is the Onsager reciprocal effect m to spin- 
transfer torque. The spin current density associated with 
spin pumping at the right interface is J|j = fi{g'j^ -|- 
(7flnx)h/47r. The annihilation of the topological charge 
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pumps spin current polarized in the z direction to the 
right metal: 

Ts _^9Rr _ f>-9R9LfJ. 

« 4 + + 

This is a central result of our work. Note that the spin 
current decays algebraically as a function of the ferro- 
magnet’s length L, which is similar to superfluid spin 
transport in an easy-plane ferromagnet [S] , but contrasts 
with the exponential decay of the spin transport by ther¬ 
mal magnons |26] . The formalism that we have developed 
is general enough to be readily extended to other easy- 
plane magnets, e.g., the case of an antiferromagnet with 
an additional easy-axis anisotropy within the easy-plane 
is closely analogous [S]. 

Brownian motion. —Let us provide an explicit model 
for Brownian motion of DWs following Ref. [^. We 
assume the following free energy for the ferromagnet: 
C/[n] = J dV{A\dxn.\^ + Kn^ — «;n^)/2, where A repre¬ 
sents the exchange stiffness, and the positive coefficients 
K and K parameterize the anisotropy magnitudes. In 
equilibrium, the local spin density s = sn lies in the xy 
plane, which can be parametrized by its azimuthal angle 
(j). A static DW solution centered at X is given by 

cos[(/>o(a: — AT)] = ±tanh[(a: — vA)/A], (6) 

with the chirality of the DW given by sgn(cos(^o • 
sin(/)o)|x>X 5 the energy Eq = 2S'/Aji, and the width 
A = \/AJk [5]. We assume here and hereafter that 
the ambient temperature is much lower than the order¬ 
ing temperature, T T^, for which thermally-induced 
changes of DW properties can be ignored. Figure [^de¬ 
picts possible types of DWs. 

The dynamics of n at a finite temperature is described 
by the stochastic Landau-Lifshitz-Gilbert (LLG) equa¬ 
tion, 

s(l-I-anx)n = n X (h-I-, (7) 

where h = —dU/dn is the effective field conjugate to 
n and is the stochastic Langevin field [18]. The 
fluctuation-dissipation theorem relates the Gilbert damp¬ 
ing constant to the correlator of the Langevin fields; 
(/i-^(r, = 2asT5{v — Y')5{t — t'). The 

Langevin equation for the overdamped dynamics of X 
can be obtained from the stochastic LLG equation by 
the collective coordinate approach [27]: 

X = 'q-^F + , (8) 

where r/ = as J dV(dx(/>o)^ = 2asS/A is the viscous coef¬ 
ficient, F = —dU/dX is the conservative force conjugate 
to X, and = —ri~^ J dV{h^^ ■ d^n) is the stochastic 
velocity [55]. The diffusion coefficient D in the corre¬ 
lator of the stochastic velocity, {v^^{X,t)v^^{X',t')) = 


2DS{X — X')d{t — t'), is related to the viscous coeffi¬ 
cient rj according to the Einstein-Smoluchowski relation: 
D = T/rj = AT/2asS (we set ks = 1). 

Nucleation and annihilation. —In the bulk of the fer¬ 
romagnet, DWs are nucleated and annihilated always in 
pairs with opposite charges [Fig. [^a)], which preserves 
the topological charge density [53]. The source and the 
drain of the topological charge, therefore, can be located 
only at the boundaries of the ferromagnet. Following 
the reaction-rate theory [55], the injection rate of DWs 
through each boundary is given by 

J±=F±(T,/r)- 7 ±(T)p±, (9) 

where is the nucleation rate, 7 '^(T) is the an¬ 

nihilation rate per unit density, and p‘^ is the density of 
q-charged DWs. 

The annihilation rate per unit density 7 ‘^(T) is the 
characteristic velocity parametrizing the escape of DWs, 
which, we expect, does not depend on the charge of 
DWs: 7‘^(T) = l{T). Interpreting the width A as the 
mean free path of DWs yields the mean thermal speed 
7(r) ~ D(T)/A. 

The nucleation rate of DWs at each interface is 
r'?(T,/r) = v{T)e:x.p[-E\y)/T], where is the en¬ 

ergy barrier for the entering of a q-charged DW and v{T) 
is the characteristic frequency that depends on details 
of the system [SO]. We take ^(T) to be independent 
of the spin accumulation y EH, similarly to the Neel- 
Brown theory for thermal switching of magnetic nanois¬ 
lands subjected to an electrical current [35] . 

Topological spin transport. —In the presence of a finite 
spin accumulation pi = fiz in the left metal, the energy 
barrier necessary for the injection of a g-charged DW is 
given by 

E‘^ = Eq + S J dx T ■ (Sc/jz) = Eq — qSgLy/4: ( 10 ) 

to linear order in pi, from which we obtain the charge- 
dependent work W‘^ = qSgLp/i done by the spin- 
transfer torque. 

When the spin accumulation is in the positive z direc¬ 
tion, ^ > 0, the entering of DWs with the positive charge 
g = 1 is favored over g = — I. The nucleation rates are 
F=*=(T, /r) = F(T)(I -I- W^/T) to linear order in pi, where 
F(T) = T{T,pi = 0) is the equilibrium nucleation rate. 
The injection rate of the topological charge through the 
left interface is given by 

= =rL{T)SW/T-^L{T)pL. ( 11 ) 

In the bulk, the topological charge current is / = —DdxP 
from the Fokker-Planck equation (H). 

At the right interface, in the absence of the nonequi¬ 
librium spin accumulation, the nucleation rate of a DW 
is independent of the charge; the topological charge does 
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not enter the ferromagnet, but only leaves it. The topo¬ 
logical charge current is, therefore, given by 


Ir = 1r[T)pr . ( 12 ) 


The conservation of the topological charge density, 
I = II = Ir, leads to the steady-state solution with 
the uniform topological charge current, 


Po 6W 

+ V-D T 


(13) 



0 1 
TjEo 


which can be recast as Eq. ([^. 7 ^ and ^r are the average 
velocity of a topological charge to cross the left and right 
interface, respectively; D/L is the average velocity of a 
topological charge traversing the ferromagnet, which can 
be seen from D/L = L/St, where 6t is the average time 
for a DW to travel the distance L. 

The spin current density by spin pumping through the 
right interface is Jfj. = Hb'r + 5 flnx)h/ 47 r. Its z com¬ 
ponent is = z • J|j = HgR(j)/4Tr to linear order in the 
bias. In the steady state with the current I of the charge 
density — 9 a;(/)/ 7 r, time evolution of the azimuthal angle (j) 
is given by 

(p/n = I. (14) 

The resultant spin current is = hpRl/A in Eq. ([^. 

Quantitative estimates. —Eor quantitative estimates, 
let us take following parameters for YIG [3133]: the spin 
angular momentum density s = 10/i/nm^, the Gilbert 
damping constant a = 10 “^, and the stiffness coefficient 
A = 5 X 10“^^J/m. Long thin YIG strips with thick¬ 
ness t = 2 nm and width iti = 50 nm are given (by 
the dipolar energy) the shape anisotropy parametrized 
by E: = 4 X 10^ J/m3 and At = 2 x 10^ J/m^ [31]. 

The algebraical decaying of topological spin transport 
manifests clearly when the ferromagnet’s length L is 
much larger than the crossover length L* = H/y, for 
which the bulk resistance dominates the interface ones, 
Rb 3> Rl,Rr. The annihilation rate is estimated as 
7 ~ D/A, which yields the crossover length L* ~ A = 
60 nm. The Boltzmann factor is exp(—E q/T) ~ 10“^ at 
room temperature T = 300 K. 

Discussion. —The topological spin transport by diffu¬ 
sion of DWs can be experimentally detected in a hy¬ 
brid structure consisting of a ferromagnetic insulator and 
two identical strong spin-orbit coupled metals, such as 
Pt|YIG|Pt (see Fig. |^, as proposed for superfluid spin 
transport [3- Given the applied electric-current density 
J£ in the left metal, the spin-current injection into the 
right metal induces the electric-current density via 
the inverse spin Hall effect, which defines the (negative) 
drag coefficient, V = —Jf^/J^. 

Figure schematically depicts the drag coefficient as 
a function of a temperature in two cases: the presence 
and absence of an easy-axis anisotropy within the easy 


FIG. 4. (color online). Schematic plot of the drag coefficient 
R = —Jr/J/ normalized by Vo{La/L) as a function of the 
ambient temperature T/Eq, where Eq is the DW energy. The 
red straight line is for the absence of anisotropy within the 
easy plane, k = 0. The blue solid curve is for the presence 
of an easy-axis anisotropy, k > 0 at low temperatures T < 
3Eo/4:. The blue dotted line shows speculative extrapolation 
of the solid curve to higher temperatures T ^ Eq. 


plane, k = 0 and k > 0. We focus on sufficiently long 
magnetic wires, for which algebraic decaying is promi¬ 
nent: L ^ Lq, for K = 0 and L ^ L* for k > 0 , where 
La = hg/2'Kas (^ 1 pva. for YIG) and g is the real part 
of the effective mixing conductance of the metal |3- For 
K = 0 , superfluid spin transport is sustained by a planar 
spiraling texture of the magnetization. The drag coef¬ 
ficient is independent of a temperature; D = Do{La/L) 
with I?o ~ 0.1 for 1 nm thick platinum [5]- For k > 0, 
superfluid-like spin transport is realized by Brownian dif¬ 
fusion of DWs. Using D = TA/2asS, the drag coefficient 
is V[T) = Tr^ApoiT)VoiLa/L) for dilute DWs, T < Eq. 
The density of DWs is given by po(F) = A~^^/KEo/ttT 
[35], which yields the blue solid line in Fig. When 
Eq ^ 0, V ^ 'Do{La/L); the algebraic decay is retained, 
provided that the temperature is well below the order¬ 
ing temperature T Tc so that the conservation of the 
topological charge is maintained [33] . 

Thermal magnons, which have been disregarded in our 
treatment, can influence the diffusive motion of DWs and 
the associated spin transfer. Thermal magnons interact 
with DWs and can affect the diffusion coefficient D at 
temperatures higher than the magnon’s energy gap m- 
This could be captured by modifying the diffusion coef¬ 
ficient D, which enters in our main result, Eq. 0. In 
addition, thermal magnons injected at the biased inter¬ 
face would exert a chirality-independent drag force on 
DWs within the spin-diffusion length [lUHISlISS]- The 
associated change in the proposed DW spin transport is 
quadratic order in the bias, and, therefore, the algebraic 
superfluid-like behavior of spin transport is not modified 
at the linear order in the bias. 

There exist other materials with exotic spin phases 
supporting localized excitations with conserved 
“charges,” e.g., monopoles in spin ices with mag- 
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netic charges [Ml- These deconfined monopoles diffuse 
by Brownian motion as experimentally demonstrated 
m. which leads us to envision that spin transport 
decaying algebraically may be implemented in a broader 
class of materials. 
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